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Bubbly turbulent Taylor-Couette (TC) flow is globally and locally studied at Reynolds 
numbers of Re = 5 x 10 5 to 2 x 10 6 with a stationary outer cylinder and a mean bubble 
diameter around 1 mm. We measure the drag reduction (DR) based on the global dimen- 
sional torque as a function of the global gas volume fraction a g i b a i over the range 0% to 
4%. We observe a moderate DR of up to 7% for Re — 5.1 x I0 5 . Significantly stronger 
DR is achieved for Re = 1.0 x 10 6 and 2.0 x 10 6 with, remarkably, more than 40% of DR 
at Re = 2.0 x 10 6 and a g i oba i = 4%. 

To shed light on the two apparently different regimes of moderate DR and strong DR, 
we investigate the local liquid flow velocity and the local bubble statistics, in particular 
the radial gas concentration profiles and the bubble size distribution, for the two different 
cases; Re = 5.1 x 10 5 in the moderate DR regime and Re = 1.0 x 10 6 in the strong DR 
regime, both at a g i bai = 3 ± 0.5% . 

In both cases the bubbles mostly accumulate close to the inner cylinder (IC). Surpris- 
ingly, the maximum local gas concentration near the IC for Re = 1.0 x 10 6 is w 2.3 times 
lower than that for Re = 5.1 x 10 5 , in spite of the stronger DR. Evidently, a higher local 
gas concentration near the inner wall does not guarantee a larger DR. 

By defining and measuring a local bubble Weber number (We) in the TC gap close to 
the IC wall, we observe that the crossover from the moderate to the strong DR regime 
occurs roughly at the crossover of We ~ 1. In the strong DR regime at Re = 1.0 x 10 6 we 
find We > 1, reaching a value of 9 (+7, -2) when approaching the inner wall, indicating 
that the bubbles increasingly deform as they draw near the inner wall. In the moderate 
DR regime at Re — 5.1 x 10 5 we find We ps 1, indicating more rigid bubbles, even though 
the mean bubble diameter is larger, namely 1.2 (+0.7, -0.1) mm, as compared to the Re 
= 1.0 x 10 6 case, where it is 0.9 (+0.6, -0.1) mm. We conclude that bubble deformability 
is a relevant mechanism behind the observed strong DR. These local results match and 
extend the conclusions from the global flow experiments as found by van den Berg et al. 
(2005) and from the numerical simulations by Lu, Fernandez & Tryggvason (2005). 

1. Introduction 

Theoretical, numerical and experimental studies on DR of a solid body moving in a 
turbulent flow have been performed for more than three decades. Drag reduction can be 
achieved by using, e.g. surfactants (Saeki et al. 2000; Drappicr et al. 2006), polymers (Virk 
1975; Berman 1978; Benzi et al. 2004; Bonn et al. 2005; White & Mungal 2008; Procaccia 
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et al. 2008) and bubbles (Madavan et al. 1984, 1985; Clark III & Deutsch 1991; Latorre 
1997; Kato et al. 1998; Takahashi et al. 2000; van den Berg et al. 2005; Murai et al. 2005; 
Deutsch et al. 2006; Sanders et al. 2006; van den Berg et al. 2007; Elbing et al. 2008; 
Gutierrez- Torres et al. 2008; Murai et al. 2008; Sugiyama et al. 2008; Jacob et al. 2010; 
Ceccio 2010). The proper implementation and understanding of DR through bubbles is 
specially relevant for naval applications, since it can lead to significant reduction of the 
fuel consumed by ships without adding substances into water. Despite all efforts done 
to understand the fundamental mechanisms behind this effect, a solid understanding 
of the DR mechanisms occurring in bubbly flows is still missing. Several mechanisms 
have been proposed as relevant, for example, bubble splitting (Meng & Uhlman 1998), 
bubble compressibility (Lo et al. 2006), bubble deformability (van den Berg et al. 2005; 
Lu et al. 2005) and effective compressibility of the flow (Ferrante & Elghobashi 2004). 
The different mechanisms as listed here, do not necessarily exclude each other and are 
likely to be active simultaneously. Depending on the Reynolds number Re - to be defined 
in §2 -, the bubble size, the bubble deformability and other parameters, one mechanism 
could dominate over the other(s) (Merkle & Deutsch 1992; Moriguchi & Kato 2002; 
Shen et al. 2006). A frequently occurring categorization on DR studies as found in the 
literature, is based on the bubble size. Studies on (sub)micron-sized bubbles, commonly 
injected in or near the boundary layer of interest, have received a lot of attention in 
the last decades (see, e.g. Madavan et al. (1984, 1985); Merkle & Deutsch (1989); Clark 
III & Deutsch (1991); Merkle & Deutsch (1992); Kato et al. (1998); Meng & Uhlman 
(1998); Kodama et al. (2000); Kanai & Miyata (2001); Moriguchi & Kato (2002); Ferrante 
& Elghobashi (2004); Murai et al. (2005); Sanders et al. (2006); Shen et al. (2006); 
Gutierrez- Torres et al. (2008) amongst many already cited references). The experimental 
results as presented in the current work uses millimeter-sized bubbles injected in the free 
stream. 

The conventional systems for studying bubbly DR are channel flows (Kodama et al. 
2000), flat plates (Merkle & Deutsch 1989; Latorre et al. 2003; Sanders et al. 2006), and 
cavity flows (see the review article by Ceccio 2010). A general explanation of the physical 
mechanism behind DR is complicated since the flow in each type of the aforementioned 
systems can exhibit cither spatial dependence like developing boundary layers, or tempo- 
ral dependence like air film build-up and bubble shedding. This and other factors make 
it difficult to collapse all bubbly DR data from different setups, see again the review 
by Ceccio (2010) on experimental work on bubbly DR measured on plates and on the 
boundary of channel flows. More recently, experiments on TC flow have been conducted, 
overcoming this problem since statistically steady states are easy to achieve (Djcridi 
et al. 2004; van den Berg et al. 2005; Murai et al. 2005; van den Berg et al. 2007; Mehel 
et al. 2007; Sugiyama et al. 2008). The TC system has the advantage of being a closed 
system and thus having a well defined energy balance (see e.g. Lathrop et al. (1992) or 
Eckhardt, Grossmann, & Lohse (2007)) in analogy with the Rayleigh-Benard system 
(Ahlers, Grossmann & Lohse 2009). In addition, TC flow offers the possibility of mea- 
suring global torque and gas volume fraction with good accuracy. Taking these various 
aspects into account, TC flow is ideal for studying the drag that a surface experiences 
when moving in a bubbly flow. 

Bubbly TC flow has been experimentally studied (Murai et al. 2005; Mehel et al. 
2007) and theoretically modeled (Sugiyama et al. 2008) for low Re numbers (< 10 4 ). The 
mechanism responsible for the DR in that regime is that the injected bubbles disrupt the 
coherent flow structures, leading to DR. The efficiency of this particular mechanism gets 
lost for high Re numbers when typical fluid velocity fluctuations are larger and coherent 
structures are increasingly unsteady. 
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In the high-Re number regime, the experimental work by van den Berg et al. (2005) 
on TC flow employing bubbles on one hand and buoyant non-deformable particles on the 
other, has proven the existence of two different DR regimes within 7 x 10 4 < Re < I0 6 . 
This work focused only on the flow's global response under influence of bubble/particle 
injection, treating the changes to the local flow parameters as a black box. The division 
between the two DR regimes roughly occurs at a Weber number We ~ 1. The Weber 
number compares the fluid inertia around the bubble surface to the surface tension. In 
the regime for Re < 6 x I0 5 and We < 1, the bubbles remain close to spherical and the 
DR is moderate, found to be around 5%. The mechanism of effective compressibility of 
the mixture, suggested in the numerical work by Ferrante & Elghobashi (2004), explains 
the DR obtained in this regime. In the regime for larger Re we have We > 1 and the DR 
is significant, reaching up to around 20%. In this regime effective compressibility becomes 
less dominant and the relevant mechanism is related to the deformability of the bubbles 
(which are much larger than the Kolmogorov length scale rjx)- In support of this view, 
van den Berg et al. (2005) find that using non-deformable buoyant particles of the same 
size leads to little DR. This interpretation has also been suggested in the numerical work 
by Lu, Fernandez & Tryggvason (2005). From their simulations these authors conclude, 
that the deformable bubbles close to the wall push high vorticity regions far away from 
the wall, reducing the drag. They suggest that this pushing is caused by the lift force on 
the bubble, which strongly depends on its shape (Magnaudet & Eames 2000). 

In a later work on the same TC setup, van den Berg et al. (2007) compare the bubbly 
DR for smooth versus rough walls within 1 x I0 5 < Re < 4 x 10 5 , again focussing on the 
flow's global response. Their results support the claim that bubbly DR is a boundary 
layer effect as bubbly DR totally vanishes for the rough-wall case. 

Significant bubbly DR with smooth walls has been measured by van den Berg et al. 
(2005) in a narrow Re regime from 6 x 10 5 to 1 x I0 6 . Given above interpretation on 
bubble deformability, it is expected that a system capable of larger Re numbers should 
lead to more DR. This assumption needs to be validated. Furthermore, previous DR 
measurements were based on the global torque, leaving many issues of the mechanism 
responsible for this effect unknown: How do bubbles modify the liquid flow? How do they 
distribute themselves inside the gap? What is the local Weber number of the bubbles? 
To gain more insight, local information on the flow and on the bubbles is required. 

In the present work we explore the Re parameter space with inner cylinder rotation up 
to Re = 2.0 x 10 6 and we measure the DR as function of the global gas volume fraction 
otgiobai- F° r Re = 5.1 x I0 5 and 1.0 x 10 6 we investigate the local modifications to the 
azimuthal velocity of the fluid due to the bubbles, we measure the local distribution of 
the bubbles and their sizes, and we present their local Weber number. These two Re 
numbers are chosen because the smaller one corresponds to the moderate DR regime 
and the larger one to the strong DR regime. 

The outline of the paper is as follows. The experimental setup and global measurement 
techniques are introduced in §2, the local measurement techniques and their accuracy in 
§3. The results on the global parameters are presented in §4, and the local parameters in 
§5. Finally, conclusions are drawn in §6, followed by further discussions and an outlook 
in §7. 
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Figure 1. Schematic view of TC flow and its important parameters. In this work we only focus 
on the torque acting on the middle section of the inner cylinder (IC), shown in blue. The outer 
cylinder (OC) is kept at rest. 



Ti T L L m id F Tj Vgap 

(m) (m) (m) (m) (m 3 ) 
0.2000 0.2794 0.927 0.536 11.68 0.716 0.111 

Table 1. Geometric parameters of the present TC system; the T 3 C. Here, ri and r are the 
radii of the inner and outer cylinders, respectively, L is the length of the full IC, L mi d is the 
length of the middle section of the IC and V gap is the total volume of the TC gap. The aspect 
ratio is given by T = Lj(r — ri) and the radius ratio by rj = Tijr . 



2. Experimental setup and global measurement techniques 

2.1. Experimental setup 

The experiments are performed in the Twente turbulent Taylor-Couette (T 3 C) facility, 
see van Gils et al. (2011a) for a detailed description of the system. Here we only give a 
brief introduction on the system parameters. We refer to figure 1 for a schematic view 
of the setup and to table 1 for the values of the geometric parameters. The IC of the 
T 3 C is divided into three separate sections. The height of the middle section of the IC 
is L m id — 0.536 m, being 58% of the full cylinder height L = 0.927 m. The IC has a 
radius of = 0.2000 m and the outer cylinder (OC) a radius of r Q = 0.2794 m. The 
maximal inner and outer angular velocities are uii/2i: = 20 Hz and uj /2tt = ±10 Hz, 
respectively. We use the radial distance r, that transforms into the non-dimensional gap 
distance f = (r — ri)/(r a — ri) with indicating the IC wall and 1 the OC wall. The 
system is fully temperature controlled through cooling of the upper and lower plates. In 
the present work we only rotate the IC and we use decalcified water as the working fluid 
and filtered instrument air for the bubble injection. 

The Reynolds number Re for pure IC rotation flow is given by 



Re = U t (r - r^)jv 



(2.1) 
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Figure 2. Schematic view of the setup showing the overflow channels through which excess 
water and air exit. The ejected water is collected in the bucket, that is continuously weighed by 
the balance, providing a g i b a i- 

where Ui — Lo^i is the azimuthal velocity of the IC wall, u>i is its angular velocity and v 
is the kinematic viscosity of the liquid. 

2.2. Torque measurement 

As shown in figure 1, the IC of the T 3 C is divided into three separate sections (top, middle 
and bottom) , each designed to measure the torque by virtue of a load cell utilizing strain- 
gauge deformation. Each section is basically a hollow cylinder suspended onto the IC's 
drive axle by two low-friction ball bearings. A metal arm connects the drive axle to the 
inner wall of each IC section. Each arm has a split in the middle that is bridged by a 
load cell (see figure 7 in van Gils et al. 2011a). The torque measurements are performed 
on the middle section of the IC, that were shown to have minimal end plate effects for 
the turbulent case of pure IC rotation (van Gils et al. 2012). 

2.3. Bubble injection and global gas volume fraction measurement 

Eight upright bubble injectors are equally distributed around the perimeter of the bottom 
plate. Each injector consists of a capillary, housed inside a plug ending flush with the 
interior wall. In all our measurements we employed the same capillaries of 1 mm in 
diameter. The gas injection rate is regulated by two mass flow controllers, resulting in 
216 liters per minute at the maximum injection rate. The injected bubble size depends 
only on the shear stress in the flow for the Re range examined here (Risso & Fabre 1998). 

When air injection is activated, the excess of liquid and air needs to go out of the 
tank to avoid pressure build up. Water and air escape through four exits - referred to 
as overflow channels - connecting the inner volume of the TC, denoted by V gapi to the 
outside. The overflow channels are located on the top plate, close to the IC. Hoses are 
attached to them and lead to a bucket, where the escaping water is collected. A balance 
registers the weight of the suspended bucket over time at a sampling rate of 5 Hz, see 
figures 2 and 5. 
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We set the weight of the bucket for single phase flow as the weight zero-offset. Once air 
injection is activated and the two-phase flow achieves stability over time, i.e. the bucket's 
weight is constant over time, the air volume contained inside V gap can be calculated. 
Considering the continuously acquired temperature of the flow, the density of the water 
at each moment is known and thus the gas volume in the tank. The global gas volume 
fraction a g i b a i is given by the gas volume in the gap as percentage of V gap . 

The top, middle and bottom 1C sections have small gaps in between them to decouple 
the torque acquisition of the middle section to the neighbouring sections. These small 
gaps introduce an increasing error over time in the global gas volume fraction estimation. 
The centripetal force pushes the bubbles through the separations and into the cavities 
between the sections. Part of the water that initially was inside, is pushed eventually 
into the flow and contributes to the fluid collected in the bucket, thus overestimating 
the actual global gas volume fraction. We conduct experiments to quantify a maximal 
error of 0.5% for a global gas concentration of 3%, and the error is smaller for smaller 
gas concentrations. 



3. Local measurement techniques 

3.f. Azimuthal liquid velocity measurement - LDA 

Laser Doppler Anemometry (LDA) is used to study the fluid's azimuthal velocity com- 
ponent of the single and two-phase flows. We make use of a LDA system from Dantec 
Dynamics, consisting of a back-scatter probe head (85 mm probe 60X81) with a beam ex- 
pander (55X21) of 1.98 expansion and an achromatic front lens of 500 mm focal distance. 
This results in a measurement volume of dimensions width x height x depth = 0.07 mm 
x 0.07 mm x 0.3 mm. The probe head is mounted on a custom-built computer controlled 
traversing system, rigidly attached to the T 3 C frame. The LDA system further consists 
of two photomultiplicrs (57X18) and a burst spectrum analyzer (BSA F80). The seeding 
particles we use are round and have a density of 1.03 g cm -3 and a mean diameter of 5 
/jm (Polyamid seeding particles, PSP-5, Dantec Dynamics). As shown in van Gils et al. 
(2012), a force balance estimation tells us that centrifugal forces are at least 4 orders 
of magnitude smaller than the drag forces acting on the particles, and hence centrifugal 
effects are negligible on the particles in our examined flows. 

When applying LDA to our setup, the two laser beams in the horizontal plane refract 
in a non-linear way due to the curved surface of the OC. Therefore the angle between 
them, 9, depends on the radial position of the focal point when located inside the gap 
between the cylinders. To correct for the curvature effect, not taken into account in the 
Dantec software, the measured speed is multiplied by the factor 

n air sin(6> air /2) 
6 = ^Tfl 7oS ' > 

^watcrSm^f/ wa ter / ^) 

where n a ; r and n water are the refractive indices of air and water, and 6* a i r and 6> watcr 
are the angles between the beams when found in air or water. We refer to the work of 
Huisman et al. (2012 b) for a more detailed explanation of the LDA method used in our 
TC setup. 

3.2. Bubble statistics measurement - Optical fiber probe 
3.2.1. Introduction 

To retrieve local information on the bubble statistics inside the flow, such as the local 
gas concentration and horizontal bubble chord lengths, we employ the optical fiber probe 
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Figure 3. (a) Microscopic photograph of the fiber tip. (b) T 3 C gap - top view: mounted 
optical fiber probe to measure the bubble statistics close to the IC wall. 



technique. This technique has been used widely in bubbly flows (Mudde & Saito 200f ; 
Xue et al. 2003; Guet et al. 2003; Luther et al. 2004; Wu et al. 2008; Xue et al. 2008; 
Martinez Mercado et al. 2010; van den Berg et al. 2011). The optical probe is a local 
phase detection device, i.e. able to discriminate between water and air in our case. Light 
that is coupled into one end of the optical fiber, will be internally reflected back at the 
other extreme; the "tip". The reflection coefficient depends primarily on the refractive 
index of the phase in which the fiber tip is submerged. By acquiring the amount of 
reflected light over time, one can obtain the bubble statistics. For a detailed description 
of the working principle, we refer to Cartellier (1990); Cartellier & Achard (1991). We 
also note that optical probes in the flow obviously disturb the flow characteristics and 
statistics. For the high Reynolds numbers here and the small probes we consider this 
problem to be minor. For the different flow geometry of the Twente water channel (but 
similar Reynolds numbers) Rensen, Luther <fc Lohse (2005) did not find any effect of 
the optical probes on the spectra, but clearly the issue of the probe back-reaction on the 
flow requires attention, in particular for closed flow geometries. 

For our measurements we use a custom made glass fiber probe (0.37 NA Hard Polymer 
Clad Multimode Fiber, Thorlabs Inc., diameter of 200/im), whose tip is stripped off the 
surrounding cladding and is heat treated to have a smooth U-shaped dome, see figure 3(a). 
The internally reflected light is collected by a photodiode and is electronically amplified 
(4-point probe electronics, Kramer's Laboratorium, Delft University of Technology). In 
all our experiments an acquisition card (National Instruments 6221) samples the fiber 
signal voltage with 16 bit resolution at a sampling rate of 120 kHz. The choice for this 
specific sampling rate was determined by precursory tests and ensures that the bubble- 
probe interactions are sufficiently resolved (see figure 4) . The optical probe sticks through 
a sensor hole of the OC with the fiber tip placed inside the TC gap, parallel to the 
azimuthal axis. It is mounted on a manual traversing system with a total radial traveling 
distance of 27 ±0.02 mm along the gap, which is around one third of the gap between the 
cylinders. Figure 3(b) shows a top view of the gap and the optical probe, made possible 
due to transparent windows on the end plates of the T 3 C. 

The accuracy of this technique has been studied in water-air flows by Julia et al. 
(2005) with bubbles of diameter 2.8-5.2 mm and velocities of 0.22-0.28 m s _1 , the latter 
one order of magnitude smaller than in our experiments. The residence time is the time 
the fiber tip is immersed in one bubble and we refer to it as T^bUe- It depends on the 
bubble size, the bubble velocity, and also on the frontal location where the probe tip 
pierces the bubble. Julia et al. (2005) report a small underestimation of T^^e, that - 
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according to these authors - is caused by the local deformation of the bubble due to 
probe-induced liquid pressure over the bubble or due to being hit by the probe in itself. 
The bubble can possibly also be decelerated or deviated from its path by an upstream 
high-pressure region induced by the probe's presence. In our experiments the bubble 
velocities are one order of magnitude larger than in Julia et al. (2005), and hence they 
have a larger added mass . While this supports the notion that it is less probable that a 
bubble is decelerated or deviated from its trajectory by the probe's presence in our flows, 
proofing this claim is hard. A force balance estimation on these effects remains elusive, 
due to the complex nature of the coupled bubble-flow interactions. Also, fully spatially 
and temporally resolved flow fields around the bubbles in the Re regime examined here, 
are hard, if not impossible with the current state of experimental techniques. Hence we 
can not guarantee that the bubbles won't be deformed, deflected or decelerated due to 
the presence of the probe, and we will take this into account by amply increasing the 
error bars of our local bubble statistics measurements, discussed in §3.2.3. As will turn 
out, even with these error bars we can make qualitative conclusions on the bubble-flow 
dynamics. 

Cartellicr (1990) introduces the concept of a latency length L for a phase-detecting 
optical fiber probe, resembling the spatial resolution of the interface detection by the 
given probe. In order to appreciate the concept, we first take a look at an example of a 
bubble-probe interaction signal as acquired with our probe. Figure 4(a) shows two bubble- 
probe events taken with our T 3 C facility, operated at Re = 5.1 x 10 5 and a g i b a i = 3%. 
The signal is detrended such that volts corresponds to the fiber tip submerged in water. 
We focus our attention on the first bubble-probe interaction, starting around time t = 
0.5 ms: A clear plateau is reached in the voltage signal, indicating that the probe has 
pierced the bubble and has had sufficient residence time in the gas-phase for the tip to 
dewet completely. The time it takes for the signal to go from 10% to 90% of the mean 
voltage making up the plateau, is an experimental determination of the duration of the 
signal rise T m , as put forward by Cartellier (1990). The 10% and 90% slice levels are 
indicated by the dashed horizontal green lines in figure 4(a), with the upwards pointing 
triangles and corresponding dotted vertical lines signalling the starting and ending time 
markers of this specific T m . In a similar fashion one can arrive at an estimate for the 
signal fall time Td - commonly smaller than T m (Cartellier 1990) - indicated by the solid 
downwards pointing triangles. By taking V$ as the velocity of the air-water interface 
hitting the probe, viz. the bubble velocity, Cartellier (1990) arrives at L = T m Vi as the 
latency length of the probe. In the next section §3.2.2 we will provide an estimate of L 
for our current probe. 

The bottom line is that L (or likewise, T m and Tj) should be short enough when 
compared to the size (residence time) of the bubbles, in order to easily resolve the 
bubble-probe interactions. If the response time T m becomes too large with respect to 
the bubble-probe interaction time, either due to too small bubbles or due to too high 
bubble impact velocities, the fiber tip will not have enough time to dewet completely. 
This then leads to signals as similar to the second bubble-probe interaction of figure 4(a) 
and the full figure 4(b), in which we show bubble-probe signals of our most challenging 
locally investigated bubbly flow at Re = 1.0 x 10 6 and a g i b a i — 3% at a radial location 
near the IC, r — 0.008. We emphasize that such signals are still usable if appropriate 
signal post-processing settings are chosen. To evaluate the optical probe against another 
independent measurement technique, we will compare the estimated horizontal bubble 
diameter as obtained from the optical probe with image analysis results of high-speed 
image recordings, presented in §5.2.2. Nevertheless, an evaluation on the latency length, 
the signal post-processing settings and the resulting accuracy is in place. 
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Figure 4. Typical detrended voltage signals (red lines with dots) of bubble-probe interactions 
as acquired by our fiber probe placed inside the T 3 C facility, (a) operated at Re = 5.f x 10 5 and 
otgiobai = 3%, displaying two distinct bubbles, (b) operated at Re = f.O x fO 6 and a g iobai = 3%, 
displaying five distinct bubbles. A positive signal slope indicates dewetting and a negative slope 
indicates rewetting of the tip. By applying a threshold algorithm on the detrended signal, based 
on the signal amplitude and its time-derivative, a bubble mask is obtained (solid blue lines). 
The width of the bubble mask is taken as the residence time Tbubbie per specific bubble. The 
dashed vertical and horizontal lines and open triangles in panel (a) aid the discussion on the 
latency length L in §3.2. Note that the plotted time range is smaller in panel (b). The signal is 
acquired at a sampling rate of 120 kHz. 



3.2.2. Validation and post-processing 

We obtain the latency length L of our optical fiber probe by evaluating the acquired 
bubble signals, as will be presented in the upcoming sections. Only bubble-probe inter- 
action signals that exhibit a clear plateau as shown in figure 4(a, first interaction) are 
selected for the evaluation. After detrending the signal, we calculate the 10% and 90% 
slice levels of each of these events in order to retrieve the average signal rise time (T m ) 
with standard deviation er m , and the average signal fall time (T^) with standard devia- 
tion ad, taken over N such events at a radial distance of f ~ 0.1. Finally, we look ahead 
at the time-averaged measured azimuthal liquid velocity of the two two-phase flow cases 
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N (T m ) Cm (Td) Od Vi L 

(ms) (ms) (ms) (ms) (m s~ ) (mm) 
Re = 5.1 x 10 5 , a g i obal = 3 ± 0.5%: 685 0.31 0.09 0.10 0.01 2.9 0.9 
Re = 1.0 x 10 6 , a global = 3 ± 0.5%: 354 0.25 0.08 0.10 0.01 5.7 1.4 

Table 2. Bubble detection characteristics of our optical probe for the two two-phase flow cases 
as studied locally in the current work. A selection of N bubble-probe interactions that lead 
to a complete dewetting of the tip is used to estimate the average signal rise time (T m ) with 
standard deviation a m and the average signal fall time (Td) with standard deviation ad- The 
bubble- water interface velocity is Vi and L — T m Vi is the latency length of the optical probe 
(Cartellier 1990). 



at that radial location (see figure 8), and we set V* equal to this velocity by assuming 
no slip or drift of the bubbles. This assumption will be examined in future work. The 
results are summarized in table 2. For the case of Rc = 1.0 x 10 6 the bubble velocities 
are larger, leading to less occurrences of complete dewetting of the tip and leading to 
smaller T m values. The signal fall time Td, however, is unaffected by the higher bubble 
impact velocities for the cases examined here. We can only speculate on the cause for 
this. The resulting latency lengths L are different for the two Re cases and are of order 1 
mm, similar to figure 12 of Cartellier (1990), which shows the same order of magnitude 
and trend for the U-shaped probe with larger latency lengths at higher Vi, although the 
interface velocities examined there are two to five times smaller than in our case. 

Evidently, the latency length of our probe is of the same order as our bubble size. 
This causes many wedge-shaped bubble-probe interaction signals with non-constant local 
signal maxima and minima, like as in figure 4(b). In this figure the first four bubbles are 
in such short succession after each other, that the tip does not rewet completely and, 
consequently, that the signal does not drop to zero before a new bubble is being pierced, 
resulting in a 'jagged' signal. Clearly, simple signal amplitude threshold schemes, as 
discussed in detail by Cartellier (1990), are not able to resolve the signal into separate 
bubble events. The post-processing algorithm we apply, makes use of the first time- 
derivative of the signal to discriminate between tip dewetting and tip rewetting. The 
complete algorithm is as follows: first, a single signal amplitude threshold is applied on 
the detrended signal in order to detect a first signal rise above the noise level. This point 
in time is set as the first bubble arrival time Tat- The noise level in all our runs was ±0.04 
V and the applied signal amplitude threshold was Vthres = 0.08 V, set as low as possible 
while still safely away from the noise level. Second, we calculate the time-derivative of 
the signal and find the first point in time T\ after Tat at which the signal rate drops 
below —8.4 x 10 3 V s -1 . This threshold value is at 75% of the signal rate associated 
with the air-film drainage as obtained from the Td study. The 75% is still large enough to 
neglect signal oscillations that sometimes happen in the case of 'overshoot' when reaching 
a signal plateau, as explained by Cartellier (1990). We also calculate the first point in 
time T2 after Tat at which the signal amplitude drops below Vthres again. The bubble 
exit time is now set to Tgx = uiin(T 1; T 2 ). Third, we set the next bubble arrival time 
to that point in time at which the signal has a positive slope combined with a signal 
amplitude larger than Vthres- The algorithm repeats itself from of point two, until the 
complete time series is processed. The bubble residence times per bubble event are now 
given by T^^Ue — Tet — Tat, identical to the widths of the corresponding binary bubble 
mask, shown as the solid blue lines of figure 4. Manual visual inspection of the algorithm 
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shows a correct single-bubblc-detection success rate of > 95%. The remaining 5% are 
bubble events that perhaps should have been split in multiple separate bubble events, or 
vice versa. 

3.2.3. Accuracy 

What can be said on the accuracy of the measurement technique? Obviously, there are 
many possible errors that are hard to quantify, but we can capture the most prominent 
ones. 

First of all, we must correct for the chance distribution of the frontal piercing location, 
leading to detected bubble chord lengths smaller than the major bubble chord. While 
this will not influence the detected local gas volume fraction, it will underestimate the 
bubble sizes. We quantify this underestimation in a first order approximation by assuming 
a spherical bubble of radius R. The mean chord length is then given by 4/3-R, which is 
1.5 times smaller than the diameter of the sphere 2R. Hence, we can expect a similar 
contribution to the underestimation of detected bubble sizes by a factor of £, c hord ~ 1-5, 
as long as the bubbles can be sufficiently approximated by spheres (see figure 7 for 
photos of the studied bubbly flows). We will take £, c hord into account when calculating 
and presenting bubble sizes. Studies on the (time-resolved) shape of individual bubbles 
fall outside of the scope of the current paper. 

We can distinguish between two different error sources. One error source comes from 
the masking algorithm. Its ability to correctly detect single bubble events in time is key 
here. With only 5% of wrongly segmented Tbubbie times, supposedly normally distributed 
around their correct values, the statistical mode of Tbubbie as used in our upcoming- 
calculations, will change insignificantly when statistical convergence is achieved. Another 
factor concerning the algorithm is how the chosen threshold values are influencing the 
location of the retrieved time-markers Tat and Tet- This imposes an uncertainty error 
£mask on the Tbubbie estimation. To our advantage here are the steep signal rates at the 
onset of tip dewetting (signal rise) and tip rewetting (signal fall). We estimate an error 
in time of ±2 samples, equivalent to £ mas k = ±16.7 /is. The quality of the bubble mask 
can be judged in figure 4. 

Another error source originates from the invasive nature of the optical fiber probe as 
studied by Julia et al. (2005), possibly deflecting, decelerating and/or deforming the bub- 
bles. As stated before, we are inclined to view the deflection and deceleration errors as 
negligible due to the large added mass of the bubbles. Nevertheless, we set the underesti- 
mation error in the local gas volume fraction due to bubble deflection to Edefiect = +5%. 
Julia et al. (2005) report an overall overestimation of Tbubbie caused by deceleration and 
an overall underestimation caused by the deformation, see figure 15 of Julia et al. (2005). 
Hence we set the error in Tbubbie due to deceleration at an overestimation of Edecei. = —5%. 
The error due to deformation we amply set at an underestimation of Sdeform — +50%. 

In the upcoming plots we will present error bars which are based on above error 
estimates. These error bars predominantly represent systematic errors. As proof for this 
statement, we refer the reader to appendix B where we present excellent repeatability 
for the local measurements at Re = 1.0 x 10 6 . 

4. Global DR measurements 

4.1. Measurement procedure 

We perform measurements on the global torque r on the middle section of the IC as 
function of the global gas volume fraction a g i b a i with the IC rotating at a fixed rate. We 
carry out experiments at Re = 5.1 x 10 5 , 1.0 x 10 6 and 2.0 x 10 6 , corresponding to IC 
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rotation frequencies of u>i /2tt = 5.0, 10, and 20 Hz, respectively. We start each experiment 
with single phase flow and once temperature stability at 21 ±0.1° C is achieved, the 
acquisition of the torque is initiated and the injection of air follows. We start with a small 
injection rate ~ 1 1/min. Initially, nearly all the air bubbles become trapped in the strong 
turbulent fluctuations of the TC flow with little air escaping into the expansion vessel 
and hence the global gas volume fraction rapidly increases. Eventually, the flow becomes 
saturated with bubbles and a dynamic equilibrium settles in between the volume of air 
entering and leaving the flow, stabilizing ctgiobal- From now on we increase the injection 
rate gradually over time, forcing more air into the flow until we reach a g i b a i ~ 4%. This 
ramp on a g i b a i takes between 25 and 160 minutes, depending on the investigated Re 
number and the gas injection rate. Higher Re number flow accumulates bubbles faster 
and hence a g i b a i increases more rapidly. See figure 5 for a typical time evolution with 
the dotted vertical line representing the onset of dynamic equilibrium. 

4.2. Results 

In figure 6 the dimensional torque data, normalized by the torque of the single phase 
flow t(0), are plotted as function of a g i a b a i- Only after dynamic equilibrium has been 
reached in a g i b a i do we represent the data by symbols instead of by dashed curves. The 
general trend of each experiment, carried out at a fixed rotation rate, shows a decrease 
in dimensional torque for increasing global gas volume fraction a g i b a i , when at dynamic 
equilibrium. Furthermore, the larger the Re number, the more effective the reduction 
of the drag is. It is remarkable that there is a major DR of 40% at Re = 2.0 x 10 6 
with ctgiobai = 4 ± 0.5%, which by far was not attainable by previous measurements in 
turbulent TC flows. 

In the figure unpublished data from van den Berg et al. (2007) have been included, 
carried out on the TC facility of the University of Maryland, or UMD TC for short. 
Their data, taken at Re = 7.9 x 10 5 , lie well between the present data of Re = 5.1 x 
10 5 and 1.0 x 10 6 . Thus good agreement for bubbly DR has been found between the 
present measurements and those by van den Berg et al. (2007) . Also note the very good 
repeatability as demonstrated by the duplicated measurement series taken at Re = 5.1 x 
10 5 and 1.0 x 10 6 . 

One could argue that the reduction of torque can be explained by just a change of 
effective density and viscosity of the liquid due to bubble injection. In Appendix A we 
introduce a non-dimensional representation of the torque, trying to compensate for this 
density and viscosity dependence on the gas volume fraction. We demonstrate that this 
can not account for the DR as observed in the strong DR regime by at least a factor of 
2 in the Re = 5.1 x 10 5 case and a factor of 4 in the 1.0 x 10 6 case. The bubbles must 
modify the flow actively in addition to simply changing the fluid properties. 

A straightforward way to define the drag reduction is DR = 1 — T(a s ; o b a ;)/r(0), where 
T(a g iobai) and r(0) are the dimensional torque values of the two-phase and single-phase 
case, respectively. Using this DR definition, the net, practical effect on the torque is 
quantified. The inset of the figure shows this DR percentage for the data points at 
constant a g \ b a i = 3% and 4% as a function of Re. Van den Berg et al. (2005) find a 
pronounced DR for Re > 6 x 10 5 and a weak DR for Re less than this. The present 
measurements show a similar result: strong DR for Re = 1.0 x 10 6 and 2.0 x 10 6 and 
moderate DR for Re = 5.1 x 10 5 . We will focus on the following two cases at a g i b a i = 
3±0.5% for the local measurements: Re = 5.1 x 10 5 with DR 6% corresponding to the 
moderate DR regime and Re = 1.0 x 10 6 with DR f» 18% corresponding to the strong 
DR regime. 
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Figure 5. Typical time evolution of the global gas concentration used in the global torque 
studies, shown here for the case of Re = 5.1 x 10 5 . The sudden steep rise from a s iobai ~ 0% to 
1% reflects the accumulation of bubbles into the flow below dynamic equilibrium, i.e. the flow is 
still able to draw in more bubbles than can escape. When the dynamic equilibrium is reached, 
here after approximately 4 minutes as indicated by the dotted vertical line, the gas injection 
rate is manually increased in a quasi-stationary way. 
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Figure 6. Normalized dimensional torque r as function of a g i b a i- Closed symbols are data sets 
from the present experiments for various fixed Re numbers. The dashed lines at the start of 
these sets indicate that the bubbly flow was not yet at dynamic equilibrium. The horizontal 
bars indicate the error of the gas concentration for the present measurements. Data are included 
from van den Berg et al. (2007) carried out on the University of Maryland (UMD) TC setup 
(fully dashed curve). Inset: drag reduction defined as the dimensional torque reduction ratio DR 
= 1 - T (a global) / r (0) for various Re numbers at constant a 9 ; j, a ; = 3% (lower branch) and 4% 
(upper branch). Colored symbols are the present results, and the open black symbols are results 
of van den Berg et al. (2007) . 



5. Local measurements 

We will study the local flow profiles and bubble statistics inside the TC gap for two 
different cases, selected out of the global DR results. They are Re = 5.1 x 10 5 , falling 
into the moderate DR regime, and Re = 1.0 x 10 6 , falling into the strong DR regime, 
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(a) 



(b) 




Figure 7. High-speed images of the bubbly turbulent flow, taken through the transparent OC 
of the T 3 C facility , for the two locally examined DR regimes at (a) Re = 5.1 x 10 5 and (6) Re 
= 1.0 x 10 6 , both at a g i bai = 3 ± 0.5%. Movies are included as accompanying material. 

both at a global gas volume fraction of a g i b a i = 3 ± 0.5%. To provide the reader with a 
qualitative sense of how such flows look like, we present high-speed image snapshots of 
the two examined Re cases in figure 7. 

We will present local profiles on the azimuthal liquid velocity, estimated bubble size, 
gas concentration and Weber number, all scanned in the radial direction at middle height 
of the TC cell. Profiles at two other axial positions for the Re = 1.0 x 10 6 case were also 
taken, which we present separately in Appendix B. 



The global gas volume fraction is maintained at a constant a g i b a i = 3.0 ± 0.5% for both 
cases of constant Re = 5.1 x 10 5 and 1.0 x 10 6 , corresponding to IC rotation rates of 5.0 
and 10 Hz, respectively. The temperature of the liquid is maintained at 21 ± 0.5°C. 

The local gas concentration profiles are scanned with the optical fiber probe for 300 
s at each radial position. The measurement time of 300 s is chosen to ensure statistical 
convergence. The first radial position is as far from the IC wall as the traversing system 
allows, which is w 27 mm from the IC wall and roughly one-third into the gap towards 
the OC. In around 20 evenly radially distributed steps we scan up to the last radial 
position, which is as close as the probe's tip can safely be to the IC wall, i.e. down to 
0.03 mm. The profile of Re = 1.0 x 10 6 is measured twice to demonstrate repeatability. 

The local azimuthal liquid velocity profiles are scanned by laser-Doppler anemometry 
(LDA) in the single and two-phase cases, and once more by particle imaging velocimetry 
(PIV) for the pure single-phase case (details to be found in Huisman et al. (2012a)). 
Both techniques are set up to provide the full radial profiles, apart from the near-wall 
boundary layers. The profiles acquired with LDA consist of around 40 equally distributed 
radial positions. The number of detected Doppler burst for each radial position ranges 
from 6 x 10 3 to 50 x 10 3 for Re = 5.1 x 10 5 , and from 28 x 10 3 to 50 x 10 3 for Re 
= 1.0 x 10 6 . The inner radial positions are increasingly harder to probe with LDA due to 
the visual shielding of the bubbles. Hence we achieve an average data acquisition rate of 
w 20 Hz close to the IC, going up to > 500 Hz when entering the bulk flow region. The 
measurement time per radial position is checked to have enabled statistical convergence. 



5.1. Measurement procedure 
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Figure 8. The radial profiles of (a) the azimuthal liquid mean velocity normalized by the velocity 
of the inner wall Ui and (b) the standard deviation of the normalized velocity fluctuations, both 
at middle height of the TC cell. The solid lines are obtained by PIV and the colored symbols 
are obtained by LDA. 



5.2. Discussion of results 
5.2.1. Azimuthal liquid velocity profiles 

Firstly, we study how the liquid velocity is modified by the injected bubbles in the two 
different Re cases by comparing the single-phase case to the two-phase case. We measure 
the time-averaged local azimuthal velocity (Ug(r)) t and the standard deviation of the 
local azimuthal velocity fluctuations u' g (r) . The results are displayed in figure 8 with (a) 
the mean velocity normalized by the IC wall velocity, which is Ui — 6.28 m/s for Re 
= 5.1 x 10 5 and Ui — 12.6 m/s for Re = 1.0 x 10 6 , and (b) the standard deviation of the 
normalized velocity fluctuations. The solid lines are the single-phase cases as measured 
by PIV and the colored symbols are data on the two-phase cases obtained by LDA. LDA 
measurements were also performed on the single-phase cases and they are found to be in 
good agreement with the PIV data: The mean within ±0.6% and the standard deviation 
within ±5% as indicated by the error bars in the lower left corner of the panels. To 
prevent the figure from cluttering up, we only show the PIV data for the single-phase 
cases. The single-phase LDA data of Re = 1.0 x 10 6 can be found in appendix B. 

The rescaled mean velocity profiles in figure 8(a) fall on top of each other in the single- 
phase flow case for both Re numbers. The bubbles clearly destroy this similarity. Since 
we do not have enough spatial resolution to resolve the boundary layers, we will focus on 
the bulk regime, taken as 0.2 < r < 0.8. The mean azimuthal velocity profiles for both 
the single and two-phase cases follow roughly the same shape in the bulk. The reduction 
of the bulk azimuthal velocity due to the bubble injection is » 4% at Re = 5.1 x 10 5 
and w 12% at Re = 1.0 x 10 6 . This is not unexpected as this velocity reduction should 
reflect the torque reduction ratio on the IC wall. 

As shown in figure 8(b), the standard deviation of the normalized velocity fluctuations 
are nearly identical in the single-phase flow case for both Re numbers. For the bubbly 
flow we observe an increase in u' s as compared to that of the single-phase flow at radial 
locations f < 0.3, i.e. close to the IC. For the higher Re case this increment is stronger. 
The next step is try to determine whether the local increase in the velocity fluctuations 
can be linked to the spatial distribution of the bubbles. 
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Figure 9. (a) Probability density functions of the bubble residence time for Re = 1.0 x 10 6 
at three non-dimensional gap distances. The mode of Tbubbie is indicated by the dashed vertical 
lines, (b) The full radial profiles of the mode. The error bars show a strong underestimation, 
that is predominantly systematical as the repeatability for Re = 1.0 x 10 6 is good. 



(a) 



•2.0 



1.5 



J 1.0 



0.5 



■ 


Re = 


5.1 


x 10 5 


• 


Re = 


1.0 


x 10 6 




Re = 


1.0 


x 10 6 


T 


tog. 


ana. 


Re = 5.1 x 10 5 


▲ 


tag. 


ana. 


Re = 1.0 x 10° 



0.05 0.10 



0.15 _ 0.20 

r 



0.25 0.30 0.35 




Figure 10. (a) Estimated mean horizontal bubble diameter obtained from the optical probe 
signals and the LDA velocity profiles. The triangles are the mean horizontal bubble diameters 
resulting from image analysis, with the error bars indicating the standard deviation. The size 
estimation from the optical probe is systematically underestimated by a factor of « 1.5 for 
both Re cases, (b) The bubble diameter expressed in the Kolmogorov length scale ijk and (c) 
expressed in the wall unit j/o- 



5.2.2. Local bubble size 

Out of the fiber probe data we compute the histogram of the bubble residence time 
Tbubbie at each radial position. Figure 9(a) shows typical probability density functions 
(p.d.f.'s) for three non-dimensional gap distances f at Re = 1.0 x 10 6 . As can be ap- 
preciated, the mode of the p.d.f.'s (defined here as the average of the bubble residence 
time corresponding to normalized probabilities higher than 0.85) is displaced the more 
to the right, the more distant the fiber tip is located from the IC. This is to be expected 
since the azimuthal flow velocity drops with increasing radius. The full radial profiles 
of the mode, modeT b ^ bbu (r), are displayed in figure 9(6). The error bars are based on 
the errors £ mas fc, Sdecel. and Edeform as discussed in §3.2.3. The most left data point for 
Re = 5.1 x 10 5 in figure 9(6) is substantially below the rest of the profile, outside of the 
error bar range. The reason for this is yet unclear to us. 
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Re 

5.1 x 10 5 
1.0 x 10 6 



(m 2 s" 3 
1.96 
13.5 



(m) 
2.63 x 10" 
1.62 x 10" 



(N m" 2 
29.7 
102 



(m s- 1 ; 

0.172 
0.320 



(m) 
5.69 x 10" 
3.06 x 10" 



Table 3. Several flow parameters as calculated out of the globally measured torque for the two 
Re numbers in the single-phase case: The energy dissipation rate e, the Kolmogorov length scale 
t)k, the wall shear stress t w , the friction velocity u T and the wall unit j/o- 



By multiplying, for each radial position, the mean azimuthal liquid velocity of the 
bubbly flow from figure 8(a) with the p.d.f.'s mode of the bubble residence time, the mean 
horizontal bubble diameter as function of the gap distance is estimated: (Dbubbiei^t = 
(Ug(r))t m0( ^ e T bubble { r ) £,chord- Several assumptions are applied here: 1) no slip or drift 
of the bubbles with respect to the surrounding liquid (to be examined in future work), 
2) the average of the instant UgT^ubbie values is equivalent to the product of {Ug)t with 
modeT^Me and 3) the shape of the bubbles is spherical enough to let the detected chord 
length be transformed to a bubble diameter by applying the systematic £, c hord factor 
described in §3.2.3. 

Figure 10(a) presents these estimated local mean bubble diameters, with the error bars 

based On Srnaski &decel. £ deform 

originating from Tbubbie- We report an overall larger 
mean bubble diameter of 1.2 (+0.7, -0.1) mm for the Re = 5.1 x 10 5 case as compared to 
0.9 (+0.6, -0.1) mm for the Re = 1.0 x 10 6 case due to the stronger shear in the latter. 
The size estimation of the bubbles is constant over the scanned radial distance for Re = 
1.0 x 10 6 , exhibiting good repeatability (see inset figure 17(6) for an exploded view of 
the data). This is a strong indication that the method has merit and that the errors are 
predominantly of systematic nature. For Re = 5.1 x 10 5 the bubble size changes weakly 
over the radial TC gap. The reason for this is unknown, but it can be covered by a small 
random error. 

We also use high-speed image recordings to measure the bubble size optically at one 
radial position, focusing at a distance f w 0.3. For Re = 5.1 x 10 5 we report a horizontal 
size of 2.0 ± 0.6 mm (mean ± standard deviation) and a vertical size of 1.6 ± 0.5 mm, 
taken over 234 in-focus bubbles. For Re = 1.0 x 10 6 the horizontal size is 1.4 ± 0.4 mm 
and the vertical size is 1.1 ± 0.3 mm, taken over 144 in-focus bubbles. The triangles in 
figure 10(a) represent these horizontal bubble sizes with the vertical error bars as the 
standard deviation. It shows that the bubbles in the lower Re case are less homogeneous 
in size as its standard deviation is larger than that for the higher Re case. It also confirms 
the estimated horizontal bubble diameter from the optical fiber signals within the under- 
estimation error bars. For both Re cases the optical probe underestimates the horizontal 
bubble size at r w 0.3 by a factor of m 1.5. Finally, the mean aspect ratio, taken as 
the horizontal over the vertical bubble size, is 1.26 and 1.30 for Re = 5.1 x 10 5 and Re 
= 1.0 x 10 6 , respectively, the latter exhibiting stronger overall bubble deformation. We 
use the term 'overall' to stress that these aspect ratios are averaged over many different 
bubbles and, hence, do not reflect the time-evolution of the shape of a single bubble. 

We now compare in figures 10(6) and 10(c) the mean horizontal bubble diameter to 
the Kolmogorov length scale r\K and the wall unit yo, both calculated out of the globally 
measured torque of the single-phase flow at the corresponding Re numbers. Table 3 lists 
these values, together with other relevant flow parameters given for completeness, see 
Lathrop et at (1992); Eckhardt, Grossmann, & Lohse (2007) for the definitions. In the 
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radial region of 0.05 < f < 0.35 the horizontal bubble diameter in the moderate DR case 
is around 46 (+25, -4) rjx and 211 (+114, -19) y®, and in the strong DR case it is around 
56 (+34, -9) ijk and 298 (+180, -46) yo, in line to what has been reported by van den 
Berg et al. (2005). 

5.2.3. Gas concentration profiles 

We calculate the time-averaged local gas volume fraction as a function of the radial 
position out of the bubble mask (see §3.2) by summing over all bubble residence times 
and dividing this by the total measurement time, i.e. 



with Ntubbies the total number of detected bubbles. The resulting radial profiles on the 
gas concentration are shown as percentages in figure 11(a). The error bars are based on 
Emost, £decei., £ deform and £ deflect-, resulting in a mainly systematical underestimation. 
Once again we stress the good repeatability as made apparent by the two independently 
measured and nonetheless collapsing profiles at Re = 1.0 x 10 6 . 

Figure 11(a) clearly shows that the bubbles tend to accumulate close to the IC bound- 
ary regime for both Re numbers. This phenomena is not surprising: In centrifugally 
driven flow there is, on average, a favorable pressure gradient pushing the bubbles to- 
wards the center of rotation, here with an acceleration a cent (r) = Ug jr equal to ~ 3g for 
Re = 5.1 x 10 5 , and ~ 9g for Re = 1.0 x 10 6 (g denotes the gravitational acceleration). 
For this estimate we took the measured azimuthal liquid velocity at the center of the TC 
gap. 

Next, we discuss the radial distance from the IC wall, r Qmax , at which the gas concen- 
tration is maximum, a max = max{(a(r)) t }. For Re = 5.1 x 10 5 this distance is J"a maJt ~ 1.2 
mm, equivalent to f ss 0.016. For the higher Re case it is located at r Qmax sa 0.6 mm, 
equivalent to r 0.008. The value of r aini _ x for each Re number is close to the respective 
mean bubble diameter as presented in figure 10(a). However, from our measurements we 
cannot conclude that the r amax we measure, is the correct radial distance at which the 
true maximum gas volume fraction occurs along the gap. It could be the case that the 
fiber probe, when located very near the IC wall, creates an up-flow high-pressure region 
that deflects the bubbles from their path, hence underestimating (a(r)) t . In this region 
the corresponding Sdefiect error might be of higher value than as set for the rest of the 
scanned radius. Further experiments are required to measure this aspect in detail. 

More importantly, we compare the amplitude of the maximum local gas concentration 
for these two Re numbers. For Re = 1.0 x 10 6 we find a max w 11 (+11, -3) %, whereas 
for Re = 5.1 x 10 5 it arrives to about 25 (+22, -3) %, which is 2.3 times higher. We 
evaluate the local Froude number to try to understand this difference. 

If the flow were fully laminar then all the bubbles would be pushed against the IC wall 
due to the radial pressure gradient associated with the centrifugally driven flow. With 
the addition of sufficiently strong turbulent liquid velocity fluctuations, the bubbles will 
experience pressure fluctuations, leading the bubbles to get spread also away from the 
wall into the bulk of the flow. We estimate the ratio of the pressure fluctuations' induced 
acceleration acting on the bubbles a press fi uc t to the centripetal acceleration a cen t{r) by 
introducing a centripetal Froude number, 
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Figure 11. Radial profiles at middle height of the two examined DR regimes: moderate DR 
at Re = 5.1 x 10 5 and strong DR at Re = 1.0 x 10 6 , both at a global gas volume fraction 
ctgiobai = 3 ± 0.5%. (a) Local gas concentration with a predominantly systemical underesti- 
mation, (b) Local centripetal Froude number with a corresponding predominantly systematical 
overestimation. The repeatability for Re = 1.0 x 10 6 is excellent. 



with u' the standard deviation of the liquid velocity fluctuations, Ug the mean azimuthal 
liquid velocity, Debbie the bubble diameter and r the radial position to be considered. 

For u' and Ug we take the measured azimuthal liquid velocity profiles of the two- 
phase LDA data of figure 8, u' s (r) and (Ug(r)) t , respectively. Burin et al. (2010) report 
the same order of magnitude for the three cylindrical components of the liquid velocity 
fluctuations inside turbulent single-phase TC flow as measured by multi-component LDA. 
At this moment we can only assume that this similarity can also be applied to the two 
two-phase flow cases, but it should clearly be tested (the radial and axial components are 
hard to measure accurately as the flow field is dominated by the azimuthal component). 
For D bubb i e we take the time-averaged radial profile (D bubble (r)) t as shown in figure 10. 
Because we use the local time-averaged profiles of each of these quantities, we retrieve 
the local time-averaged Froude number (Fr cen t(r))t, presented in figure 11(b). The error 
bars are based on the small random errors in (U$(r))t and u' g {r), together with e mas k, 
£decei. and Sdeform, arising from the (D bubb i e (r))t estimation. The errors in the bubble size 
estimation are dominant and hence Fr cent is mainly systematically overestimated. For 
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both Re cases we find (Fr cen t(r)) t > 1, also within the error bars. Hence, the turbulent 
fluctuations are strong enough to diffuse the bubbles also against the centripetal force 
towards the bulk region of the flow, away from the inner wall. We do not claim that this 
evaluation is exact and it should be regarded as a first approach to try to explain the 
gas distribution differences. 

Now we evaluate the local (Fr cent (r))t number at the radial position r Qmax to under- 
stand the different a max for each Re case. The closest we can get is f = 0.06 for which 
wc find (Fr oent ) t « 1.6 (+0.4, -0.7) for Re = 5.1 x 10 5 and (Fr cent ) t 3.4 (+1.1, -1.5) 
for Re = 1.0 x 10 6 . The lower Re case has a lower Fr cent number by a factor of ps 2.1, 
indicating that the effective centripetal force on the bubbles is higher and hence, the 
bubble accumulation is stronger near the IC wall at the expense of a lower concentration 
in the bulk, as observed. The Fr cent number ratio at this radial position matches the 
maximum gas concentration ratio of ~ 2.3 between both Re cases. 

The spatial broadening around the gas concentration maxima of figure 11(a) implies a 
wider bubble accumulation zone for Re = 5.1 x 10 5 as opposed to Re = 1.0 x 10 6 . Clearly, 
there must be more at play than just the balance between pressure fluctuations and 
centrifugal forces on the bubbles for explaining the concentration profiles: If this would 
not be the case then the Fr cent number would dictate the narrowest spatial distribution 
to happen for Re = 5.1 x 10 5 . 

What is unambiguously remarkable is the following: One would expect that the DR 
efficiency is correlated to the bubble concentration near the inner wall, and that a higher 
local gas volume fraction should give a larger DR. However, the present measurements 
show an opposite trend: The local gas concentration near the inner wall is higher in the 
moderate DR regime than that in strong DR regime. This clearly suggests that a higher 
local gas concentration near the inner wall does not guarantee a larger DR. What then 
is the dominant effect for the strong DR at the high Re numbers? To try to answer 
this question we will look into bubble deformability by defining the local bubble Weber 
number. 

5.2.4. Local Weber number 

Following van den Berg et al. (2005), we define the bubble Weber number as 

We ^ pu' 2 D bubbl% 

G 

with v! the standard deviation of the velocity fluctuations and a the surface tension of the 
bubble-liquid interface. Relation (5.3) is identical to what is used in van den Berg et al. 
(2005). Some typographical errors appear in that work: The Weber number definition 
lists the radius of the bubble R b while it should read the bubble diameter D bubb i e . This 
concerns their equation (4) and the Weber number definition three lines above it. Also, 
their abstract should read l R = 0.5 mm' as the value of the used bubble radius and in 
the conclusion the statement 'R b = 42r? = 300?/o' should read 'R b = 21rj — 150yo'- The 
rest of their values, including the right-hand side of equation (4) and their graphs are 
correct. 

For v! we take the azimuthal liquid velocity fluctuations u' g (r) of the two-phase LDA 
data, see figure 8(b), following the same reasoning as used for the Fr cent number es- 
timation in section 5.2.3. For D bubb i e we use the local time-averaged bubble diameter 
(D bubb i e (r))t of the optical fiber probe data, see figure 10(a). Substituting the values p 
= 998 kg m~ 3 and a = 0.0727 N m -1 , corresponding to water at 21°C, provides us now 
with the local time-averaged Weber number of the bubbles inside the TC gap, scanned 
over 1 /3 of the radius close to the IC at middle height of the TC cell. We present these 
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Figure 12. Local time-averaged Weber number profiles of the two examined DR regimes at a 
global gas volume fraction of a g i a t a i = 3±0.5%: Moderate DR at Re = 5.1 x 10 and strong DR 
at Re = 1.0 x 10 6 . The vertical' error bars indicate a predominantly systemical underestimation, 
originating from the bubble size estimation. 



(We(r))t profiles for the two examined DR regimes in figure 12. The vertical error bars 
are based on the small random error in u' g (r) together with s ma sk, £decei. and Sdeform, 
arising from the (D bubb i e (r)) t estimation. The errors in the bubble size estimation are 
dominant and hence (We(r)) t is mainly systematically underestimated. 

For Re = 1.0 x 10 6 we see (We(r)) t > 1 and that the bubbles increasingly deform 
when approaching the IC wall, where the liquid velocity fluctuations are stronger, finally 
leading to (We) t « 9 (+7, -2) closest to the inner wall. The bubble Weber number for Re 
= 5.1 x 10 5 is lower and it remains nearly constant with a value around 1, increasing to 
(We)t sa 1.4 (+1.0, -0.3) closest to the inner wall. This implies that the bubbles remain 
nearly spherical and/or retain their shape over time, because the velocity fluctuations are 
not strong enough to deform the bubbles, even though the estimated bubble size is larger 
- (D bubb ie(r))t ~ 1-2 (+0.7, -0.1) mm - than that for the higher Re case - (D bubble (r)) t rj 
0.9 (+0.6, -0.1) mm. At Re = 5.1 x 10 5 the global DR is around 6%, whereas the DR is 
about 18% for the case of Re = 1.0 x 10 6 , having bubbles of higher deformability. Our local 
and global results thus clearly demonstrate that the bubble deformability is important 
for the observed strong DR, which is consistent with the findings by van den Berg et al. 
(2005) and van den Berg et al. (2007), based on their global torque measurements, and 
by Lu, Fernandez & Tryggvason (2005), based on their numerical calculations. Indeed, it 
is this bubble deformation that is the origin of the lift force enhancement observed by Lu, 
Fernandez & Tryggvason (2005), resulting in the bubbles pushing away the liquid from 
the wall. Here we provide local experimental evidence on the important role of bubble 
deformability for strong DR in bubbly turbulent TC flow. 



6. Conclusion 

We have shown that the net DR, based on global dimensional torque measurements 
in bubbly turbulent TC flow with IC rotation up to Re = 2.0 x 10 6 and stationary OC, 
increases with increasing global gas volume fraction a g i oba i < 4% and that the efficiency 
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of the DR increases for increasing Re numbers. We find a remarkable net DR of 40% 
obtained at a g i oba i = 4 ± 0.5% for Re = 2.0 x 10 6 , which was not attainable by previous 
TC setups. The present global measurements at Re = 5.1 x 10 5 and 1.0 x 10 6 show 
good agreement with previous global measurements from van den Berg et al. (2007). We 
observe a moderate DR for Re = 5.1 x 10 5 up to 7%, and a strong DR for Re = 1.0 x 10 6 
and 2.0 x 10 6 up to 23% and 40%, respectively. 

To understand the flow differences between the moderate and strong DR regimes, we 
studied the flow locally for two Re numbers at fixed a g i oba i = 3 ± 0.5%, each falling into 
a different DR regime. The smaller one, Re = 5.1 x 10 5 , represents the moderate DR 
regime at DR = 6% and the larger one, Re — 1.0 x 10 6 , represents the strong DR regime 
at DR = 18%. Consequently, we have extended the global flow experiments by van den 
Berg et al. (2005, 2007) by providing local flow information. 

We investigated the azimuthal liquid velocity in the single and two-phase flows for 
each Re case. The amplitude of the mean azimuthal velocity in the bulk (0.2 < f < 0.8) 
for the bubbly flow decreases « 4% at Re = 5.1 x 10 5 , and w 12% at Re = 1.0 x 10 6 
as compared to that of the single-phase case. These reduction ratios match the global 
torque reduction ratios. Furthermore, we measured increased liquid velocity fluctuations 
close to the IC in both two-phase cases as compared to their equivalent single-phase case. 

For retrieving local information on the bubble statistics we made use of the optical 
fiber probe technique. We have discussed the possible error sources of this technique 
and we have demonstrated that it underestimates the bubble size, predominantly due 
to systematic errors, by a factor of rj 1.5 in our examined Re cases as compared to 
high-speed image analysis. 

The local gas concentration measurements show that the bubbles tend to accumulate 
close to the IC boundary regime for both Re numbers. This can be explained by the 
centrifugally driven flow, pushing the bubbles towards the IC wall with an acceleration 
of around 3 to 9 g. The velocity fluctuations are strong enough to diffuse the bubbles also 
towards the bulk of the flow, as expressed by the centripetal Froude number Fr cent > 
1. The radial distance from the IC wall at which we found the maximum local gas 
concentration is close to the respective bubble diameter. Surprisingly, the maximum local 
gas concentration at Re = 5.1 x 10 5 , i.e. in the moderate DR regime, is s» 2.3 times higher 
than that at Re — 1.0 x 10 6 , i.e. in the strong DR regime. We found agreement between 
the time-averaged local (Fr ceni (r)) t number and the maximum gas concentration near 
the inner wall, explaining this large concentration difference in a first approach. More 
importantly: This result unambiguously reveals that a higher local gas concentration 
near the inner wall does not guarantee a larger DR. 

By estimating the time-averaged local Weber number (We(r)) t of the bubbles, we 
concluded that in the strong DR regime at Re = 1.0 x 10 6 the larger liquid velocity 
fluctuations deform the smaller bubbles of {D bubb [ e (r)} t ~ 0.9 (+0.6, -0.1) mm, since 
(We(r)) t > 1. The local Weber number increases with decreasing distance from the inner 
wall, reaching (We)t ~ 9 (+7, -2) close to the inner wall, implying that the deformation 
of the bubbles is stronger near the IC wall. On the contrary: In the moderate DR regime 
at Re = 5.1 x 10 5 the bubbles remain close to spherical and/or retain their shape over 
time more, because (We(r)) t ~ 1 with near the inner wall (We); » 1.4 (+1.0, -0.3), even 
though the bubbles are larger in size: (D bubb i e (r)) t 1.2 (+0.7, -0.1) mm. 

We conclude that increasing the gas volume fraction has a positive influence on the 
azimuthal liquid velocity decrease, which is responsible for the observed DR. What is 
important for strong DR is that bubbles deform (We > 1) close to the boundary layer 
regime of the IC wall, leading to a large azimuthal velocity decrease and thus a strong 
DR. These local results match and extend the conclusion on the global torque results as 
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Figure 13. Typical time series of the measured local gas volume fraction, here post-processed 
with a sliding average window of 1 second in duration, showing the short-time fluctuations and 
long-time stability. The circles are the windowed (quasi-instantaneous) gas volume fraction, the 
solid vertical line is the mean and the dashed vertical lines are the standard deviation of the 
fluctuations around the mean. The time series correspond to the Re = 1.0 x 10 6 measurement 
(stars) at the radial location of maximal gas volume fraction r amax of figure 11(a). 



found by van den Berg et al. (2005) and van den Berg et al. (2007), and on the numerical 
results by Lu, Fernandez & Tryggvason (2005), both stressing the central role of bubble 
deformability. 



7. Discussion & outlook 

No periodicity has been found in the bubble signal time series at any of the examined 
radial locations. This includes spectral and auto-covariance analysis on bubble and inter- 
bubble arrival times. For completeness, we show in figure 13 a typical time series of the 
measured local gas volume fraction, here post-processed with a sliding average window 
of 1 second in duration, revealing the short-time fluctuations and long-time stability. 

For the Re — 5.1 x 10 5 case we occasionally saw bubbles coalesce into air patches of 
several centimeters in size. Sometimes, these air patches traveled along the inner wall, 
supporting the notion of a lubricating gas film. Mostly, however, the air patches seemed 
to be inside of the bulk, albeit with strong upward motion. For Re = 1.0 x 10 6 the 
observation of air patches is even less frequent, and if they exist then they are more often 
observed to be inside of the bulk and not in contact to the wall. These events correspond 
to the long right tails of the bubble residence time p.d.f.'s and their probability is 3 
orders of magnitude smaller than that of the mode. Hence, we are of the opinion that 
the occasional air patches do not influence our results noticeably. 

We note that the bubble distribution is a two-way coupled problem: Bubbles are drawn 
into a radial distribution as suggested by Fr cerat (r) and at the same time the magnitude 
of Fr cent (r) is determined by the velocity fluctuations induced by the presence of the 
bubbles. And likewise, the effect of bubble deformability also influences the velocity 
fluctuations' strength. The effect of deformability can be eliminated by injecting rigid 
neutrally buoyant particles into the flow, similar to what was used in experiments by 
van den Berg et al. (2005). One could also further study above problem by changing the 
Fr ce „t and We numbers independently of each other. This can be achieved by changing 
the liquid density p as it enters the We definition but not the Fr ce „ t definition. Note that 
changing the surface tension a will not result in independently changing Fr cerlt from We 
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as a different surface tension results in a different mean bubble diameter, which enters 
both the Fr cent and We definitions again. 

In this present work the bubble diameter is linked to the driving control parameter Re, 
because the shear strength of the flow dominates the bubble size (Risso & Fabre 1998). 
One can achieve a different scaling of the Re number to the bubble size by adding a low 
concentration of surfactant into the bubble carrier liquid. Even small concentrations of 
surfactant are known to significantly increase the effective surface tension of the bubble- 
liquid interface as it results in changing the bubble surface boundary condition from a 
slip to a no-slip condition, see e.g. the review article by Takagi & Matsumoto (2011). 
Different low concentrations of different surfactants can result in a variety of mean bubble 
sizes and deformability, and hence different We. Of course, the huge shear strength of our 
flows will still play a role. Surfactants also prevent bubble coalescence and lead to a more 
monodisperse bubble size distribution, both aiding the simplicity of the flow. We assume 
that adding small amounts of surfactants will not change the turbulent flow properties, 
which should be tested experimentally. The bubble deformability and the connection to 
the global torque DR can and should be studied further in bubbly turbulent TC flow. 

So far we have examined bubbly turbulent TC flow by only rotating the IC and keeping 
the OC at rest. The flow profile inside the TC gap becomes completely different once 
OC rotation is activated, see van Gils et al. (2012) for angular velocity profiles of such 
flows obtained in the same TC facility. In the co-rotation regime the flow will become 
less turbulent and is hence not the first choice to examine. In the case of counter-rotation 
the flow will experience more shear and a 'neutral line', i.e. the radial position where the 
mean angular velocity is zero, will exist inside the TC gap, thereby changing the overall 
radial pressure gradient induced by the centrifugal force. This could have a pronounced 
effect on the bubble distribution, possibly allowing indirect control of the radial position 
of maximum gas concentration r Qmax from the boundary layer region towards the bulk 
region. 

Another route to proceed is to work towards a one-to-one comparison between our 
experiments and the front-tracking numerical simulations by Tryggvason and coworkers 
for the TC geometry, which, at least for lower Re numbers, should be achievable. 
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Appendix A. Non-dimensional torque reduction ratio 

When injecting bubbles into a liquid at small concentrations, the effective kinematic 
viscosity and density of the two-phase fluid can be approximated, respectively, as follows: 



Pa g iobai - Pi 1 - a global), 
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v a global = V(l + ^global), 
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Figure 14. Normalized non-dimensional torque G as function of otgiobai- This figure shows an 
alternative representation of the data in figure 6. The solid green curve shows the suggested 
dependence of the density and viscosity terms of the normalized non-dimensional torque on 

^global • 

where a g i b a i is the global gas volume fraction, and p and v are the density and viscosity 
of the liquid at a g i b a i = 0. The 'Einstein' equation (A 2) (Einstein 1906) is a first order 
approximation only and assumes the dispersed phase to be non-self-interacting and to be 
of the same density as the carrier phase. Though the applicability becomes questionable 
for bubbly turbulent flows with up to 4% gas concentration, we will use (A 2). 

To estimate the magnitude of the torque reduction, purely caused by the change of 
fluid properties when injecting bubbles into the flow, we define a non-dimensional torque 
based on the two-phase viscosity and density as: 



G(ctglobal) 



2np, 



Cglobal OL g i 



V a q lobal Lmid 27rp(l - a g l bal)v 2 (l + 2 a global) 2 L mld 



(A3) 



with t as the measured torque on the middle section of the IC. We focus on the effect of 
bubbles on the torque and neglect the change of the Re number. 

The experimental results on the non-dimensional torque reduction are plotted in figure 
14, comparable to the dimensional torque presented in figure 6. The reason for the shift 
in the UMD data can be contributed to a different fluid temperature of 24 °C. The solid 
green curve represents the suggested dependence of the density and viscosity terms of 
the normalized non-dimensional torque on a g i bai- Therefore, the difference between the 
solid green curve and each of the measured curves tells us the contribution to the DR 
that is not solely caused by a change of fluid properties. In the case of Re = 5.1 x 10 5 
the difference remains very small even up to a g i b a i ~ 2.5%, but it increases significantly 
when ctgiobai > 2% for Re = 1.0 x 10 6 and 2.0 x 10 6 . Hence, the bubbles in the strong DR 
regime must modify the flow actively in addition to simply changing the fluid properties. 



Appendix B. Axial dependence at Re = 1.0 x 10 6 

We also investigate the axial dependence of the flow at Re = 1.0 x 10 6 with a g i bai = 
3 ± 0.5%. The OC provides three locations through which the optical probe can be 
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Figure 15. Snapshot of the T 3 C setup at Re = 1.0 x 10 6 and a g i bai = 3%. The axial mea- 
surement locations are tagged. Possible remnants of turbulent Taylor vortices are visible, that 
seem to coincide with the measurement locations. This is only coincidental at the moment of 
the snapshot. Neither stable vortices nor statistically stationary vortices have been observed in 
the bubbly flow. 

inserted. These are located at z mi ^ = L/2, zi, ot — z mi d — 15 cm and z top — z m id + 15 
cm, see figure 15. In this appendix we will omit showing the error bars on the bubble 
statistics as they are of predominantly systematical nature and, hence, do not aid the 
axial dependence investigation. 

We will first focus on the azimuthal liquid velocity profiles as measured by LDA at 
these locations: The normalized mean velocity (Ug(r)) t /Ui in figure 16(a) and the turbu- 
lence intensity u' g (r)/ (Ug(r)) t in figure 16(6). For the single-phase case we find no axial 
dependence over the examined range as both the mean velocity and turbulence intensity 
profiles are identical to within 0.7%. Hence we represent the single-phase case as the 
black stars. When we inject bubbles into the flow we observe a weak axial dependence on 
the mean azimuthal velocity but not on the turbulence intensity, see the colored symbols. 
The relative mean velocity amplitude difference in the bulk of the flow for the Zt op and 
Zbot profiles is ±4% when compared to the z m id profile. Can this be linked to the bubble 
distribution? 

So, we examine the corresponding gas concentration profiles, see figure 17(a). Two 
independent measurements are performed at each height to show the measurement re- 
peatability. We note an increase in the maximal value of the local gas concentration 
(a(r))t for increasing height. Between zt op and z m id and between z m id and Zbot there is 
a factor 2.8 difference. These increments in gas volume fraction can be due to two 
factors: More bubbles or bigger bubbles. 

Hence, we calculate, in the same manner as described in §5.2.2, the local bubble diam- 
eter individually for each axial position, see the inset of figure 17(b). No axial dependence 
on the bubble size is found. When we consecutively calculate the average bubble detec- 
tion rate (fbubbie( r ))t for each axial position, we indeed find a match to the maximum 
gas volume fractions, see figure 17(b). Interestingly, the local bubble distribution profiles 
show universality when normalized by their respective maxima a max , see the inset of 
figure 17(a). 

The reason for this strong axial gradient on the gas volume fraction becomes appar- 
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Figure 16. (a) The radial profiles at different axial locations of (a) the liquid azimuthal mean 
velocity normalized by the IC wall velocity Ui, (b) the turbulence intensity. The error bars in 
the lower-left corners indicate the errors on the data series. 
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Figure 17. (a) Local gas concentration profiles at three axial locations in the flow at Re 
= 1.0 x 10 6 and a g i i, a i — 3%. Inset: The local gas concentration profiles normalized by their 
respective maxima a max , hence revealing universality, (b) The local average bubble detection 
rate in number of bubbles Ntutbies per second. Inset: Corresponding local mean bubble diameter. 



ent when exploring the following force balance. We find that the bubbles experience a 
centripetal acceleration Ug jr ~ 9g in the bulk of the flow as compared to the lg of axial 
upward acceleration due to buoyancy. These factors are not well enough apart to neglect 
the influence of gravity on the gas distribution. 

We conclude that the local gas concentration depends strongly on the axial position. 
However, the liquid mean velocity profiles change only weakly and the turbulence inten- 
sity profiles fall on top of each other. This suggests that the azimuthal mean velocity 
reduction is not sensitive to the local gas concentration along the axial direction. 
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